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C^l ' We study spherically symmetric solutions to the Einstein field equations under the as- 

sumption that the space-time may possess an arbitrary number of spatial dimensions. The 
general solution of Synge is extended to describe systems of any dimension. Arbitrary di- 
mension analogues of four dimensional solutions are also presented, derived using the above 
scheme. Finally, we discuss the requirements for the existence of Birkhoff's theorems in 
space-times of arbitrary dimension with or without matter fields present. Cases are dis- 
cussed where the assumptions of the theorem are considerably weakened yet the theorem 
still holds. We also discuss where the weakening of certain conditions may cause the theorem 
to fail. 
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■ §1. Introduction 

O 

Space-times possessing dimension greater than four have been of much in- 
terest at least since the pioneering ideas of Kaluza 1 ) and Klein 2 ) 3 ). Since then there 
have been numerous theories of unification (for example superstring theory) , many of 
which require more than three spatial dimensions to be consistent. The low energy 
sector of many of these theories reduce to a multi-dimensional General Relativity 
', theory as is studied here. 

Although extra dimensions are usually thought to be compact, the extra 
dimensions may be manifest on scales which are relevant when studying cosmologi- 
cal systems near the big bang or gravitational collapse approaching the singularity. 
Also, there has lately been much interest in the possibility of large extra dimen- 
sions 4 )~ 6 ) in which higher dimensional effects may be observed at relatively low 
TeV scales. Microscopic systems above these scales may be formed from the collapse 
of large, effectively four dimensional, initial conditions. If these scenarios do indeed 
describe our universe, then gravitating systems at these scales will behave as higher 
dimensional systems and deviate from the predictions of four dimensional physics. 

Much interesting work has been done in the field of higher dimensional 
gravity (for example, see 7 ) - 14 ) and references therein. The excellent review by Mel- 
nikov 10 ) has an extensive list of references) . Usually, these studies involve specific 
matter fields or a class of metrics such as the FRW cosmological metrics. However, 
little work has been done on a reasonably general methodology which may be utilised 
in studying higher dimensional gravitation problems. Granted, a completely general 
method to solve Einstein's field equations does not exist and therefore in this note 
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we focus on spherically symmetric systems. Since some interesting studies have been 
performed on the subject of higher dimensional black holes (for examples see 16 ), 17 ) 
and references therein) our concentration here will be on non-vacuum systems. 

Finally, studies of General Relativity in arbitrary dimension will also serve 
to shed light on the theory's internal consistencies and lead to a greater understand- 
ing of the theory as a whole. This avenue has already proved fruitful in the case of 
low dimensional black holes 18 -*~ 22 \ for example. 

We therefore believe that it is of much interest to study how the presence of 
an arbitrary number of spatial dimensions affects the solutions of general relativity. 
To this end, in section 2, we generalise Synge's 23 ) method of solving the spherically 
symmetric Einstein field equations so that it may apply to D-dimensional systems. 
In section 3 we present solutions which are arbitrary dimensional counterparts to 
some important 4 dimensional solutions. Finally, in section 4, we rigorously prove a 
D-dimensional staticity or Birkhoff 's theorem and comment on situations when the 
theorem can fail. 



§2. The general solution 



The fundamental equations governing the space-time geometry may be de- 
rived from the action **) 



f(~ — 
J V 16vr 



+ £m I yfgdPx 



(2-1) 



where C m is the matter Lagrangian density and the positive integer D is assumed 
to be larger than 2. 

The action principle with suitable boundary conditions gives rise to the 
D-dimensional Einstein field equations: 



1 



R*i --Rd f i = 8vrT^ , 



(2-2) 



along with supplementary equations governing the behaviour of the matter fields. 
We wish to study solutions under the ansatz of spherical symmetry. In curvature 
coordinates this allows us to write the space-time metric as: 



ds 2 = _ e „(r,t) dt 2 + e X(r,t) ^2 + r 2 d[2 f D ^ 



where dQ 2 D2 ^ is the line element on a unit D — 2 sphere: 



D-3 



n=l 



\m=l 



(2-3) 



(2-4) 



*^ The literature on lower dimensional black holes is extensive. We apologize that we cannot 
cite all the excellent work in this area. 

**-* Conventions in this paper follow those of 15 ^ with Gd, the D-dimensional Newton's constant, 
and c = 1. Here, Greek indices take on values — > D — 1 whereas Latin indices take on values 
1 — > D — 1 (spatial indices) . 
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The corresponding hyper-surface volume is given by: 

_ 2^ D - 1 )/ 2 
^(D-2)-sphere = p ( jP-1 )~ ' 

(2-5) 

and the coordinate ranges are: 

h < t < t 2 , n < r < r 2 , < (o) , 0(i), . • • , 0(£>_ 4 ) < 7T, < (D _ 3) < 2vr. (2-6) 
The above metric yields the following field equations for (2-2): 



871-71 



87rT r r 



8irTt = - 



D-2 

"IP" 
D-2 

2r 2 

D-2 
1 

2r 

e -"( r .*) 



+ re 



-A(r. 



*)A(r,t), 



+ 



4 

e -A(r,t) 



(-D-3) (^l- e - A (^) x 
(£> - 3) (l - e- x ^ - re-^ r '%(r,t), r 

u(r, t), t X(r, t\ t - (A(r, t), t ) 2 - 2A(r, *), t)t 



(2-7a) 
(2-7b) 
(2-7c) 



2t/(r, t), r , r + (!/(r, t), r ) 2 + 2(Z) 3) (i/(r, t) - A(r, i)) >r 



-i/(r,*),rA(r,*), r + ^(D-3)(D-4) 



2(D-3)(D-4) 



• (2-7d) 



Enforcing conservation laws, T^ ;At = 0, yields the following non-trivial equa- 



tions: 



T% t + Ti r + l -\{r,t), t (r*-r;) + ir^ 

= o, 

L> - 2~ 



(A(r,*) + I /(r,t)) >r + 



2(D - 2) 



+ t* t + 



r; + -Kr,*) + AM)] )t rt 



-i/(r, t) r T*,+ T ft (A) 



0. 



(2-8a) 



(2-8b) 



The field equations must now be solved. We adopt here a similar method 
to that of Synge 23 ) which we generalise to accommodate arbitrary dimension. 

The system possesses six partial differential equations (which admits two 
differential identities) and six unknown functions: u(r,t), X(r,t) and the four rel- 
evant components of the stress-energy tensor. One may therefore either prescribe 
two of these functions or else they may be determined by other means, such as sup- 
plementary matter equations. We assume, for the moment, that T* and T r r are the 
known functions (we will mention more on this issue shortly). They may also be 
related by an equation of state. We solve for A(r, t) via (2-7a) by noting that (for 
D > 2) it may be written as: 



-167T 



„D-2 



D-2 



«=[r D -Mr,t)}. 



(2-9) 
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with w(r,i) := 1 - e~ A ( r >*). 
manipulation gives: 

167T 



Integrating this equation followed by minor algebraic 



-A(r,t) 



1 + 



/' 

■/rn 



T l t {x,t)x D - 



dx + 



/(*) 



1 



2m(r, i) 



„D-3 



(2-10) 



(£-2)H>- ,-„, 

Here /(i) is an arbitrary or free function of integration. To avoid a singularity at 
r = one must set this function to zero. However, here we shall keep it for generality. 

The metric function, u(r, t), is obtained from a linear combination of (2-7a) 
and (2-7b) as: 

D ~ 2 - X ^Hr,t) + X(r,t)] r , 



Svr [T\ - r r ] 
which, using (2-10), yields 



2r 



(2-11) 



o"(r,t) 



2m(r, t) 



,D-3 



exp <h(t) + 



167T 

D - I 



rr ' 
Jro . 



T r r (x,t)-T\(x,t) 



x 



D-3 



2m(x, t) 



x D ~ 2 dx 



(2-12) 



The function h(t) is a function of integration which may be absorbed in the definition 
of a new time coordinate via the transformation t = J exp [h(t)/2] dt. This is not 
always possible as will be discussed in section 4. 

The energy flux, T* , may now be defined by the equation (2-7c) and the 
lateral pressure, T 9 g , is defined from the conservation law (2-8a): 



0(A) 



D 



+ 



T r r + T 



r,t + 



1 . , D - 



-Kr,t) + A(r,t)]/ 
T r r -^(r,t\ r T\ 



(2-13) 



as this must be the lateral pressure if the energy density and parallel pressure are 
known. At this point, it can be shown that all equations and identities are satisfied. 

If one is interested in solutions given by specific matter fields, the number of 
unknowns versus the number of equations may be different than previously specified. 
For example, if the system to be studied respects absolute spatial isotropy (as in the 
case of a perfect fluid), then (2-13) becomes a differential equation for T r r = T e g 
which must be solved. The general solution presented above will still contain these 
as specific cases although other constraints must be met (which may be variationally 
derived and lead to a determinate system). 

Otherwise, one is free to prescribe the two functions. The most physical 
prescription involves specification of the energy density and one (parallel) pressure. 
This method is quite useful in examinations of relativistic stellar structure and col- 
lapse dynamics where one usually prescribes a reasonable energy and pressure from 
nuclear theory and studies of plasmas 24 ) . 



§3. D - dimensional counterparts to known solutions 



In this section we construct D-dimensional analogues of some four dimen- 
sional solutions. It is useful to study specific solutions since they are illustrative of 
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n = -^w^ - i-AK, (3-1) 



how various physical properties in the construction depend on space-time dimen- 
sion. The first solution we present is previously known. We briefly present it here to 
illustrate how it may quickly be derived using the method of the previous section. 

3.1. Kottler solution 

Consider the case of a Kottler (Schwarzschild-(anti) de Sitter) solution 25 ) 
in D-dimensions. In this situation, the cosmological constant is best viewed as part 

Of Ty\ 

87Tr D-2 V J t v 8jt 

M being a constant related to the effective mass of the black hole. Here we are loosely 
using the coordinates r and t to represent the domain within the black hole's event 
horizon as well as the exterior domain. The integrals in (2-10) and (2-12) may easily 
be evaluated to yield the D-dimensional counterpart to Kottler's original solution (a 
coordinate re-scaling is assumed): 

, , / 2M 2Ar 2 \ , dr 2 

ds = - 1 pt^t - r dt 1 + r- 

V r D - 3 (D-1)(D-2)J 2M _ 2Ar* \ 

\ L r D - 3 (D-1)(D-2)J 

+ r 2 df2f D _ 2) . (3-2) 

The term M/r D ~ 3 is a potential or harmonic function in a D — 1 dimensional Eu- 
clidean space which is related to the higher dimensional Newtonian potential. 

3.2. Homogeneous, incompressible star 

We next study the arbitrary dimensional analogue of the static, constant 
density star 26 ) as well as derive the corresponding maximum mass/radius relation- 
ship. The constant density sphere is of interest since it will yield the upper limit 
to the surface gravitational red shift for spherical (non black hole) systems in any 
dimension. Interesting examples of 2 + 1 and 4 + 1 dimensional stellar models may 
be found in 27 )~ 29 \ 

For the constant density homogeneous sphere we have 

—T\ := ( P n ° l OT r<a (3-3) 
I for r > a 

with po and a some positive constants. All principal pressures are equal and will be 
denoted by p. 

The only non-trivial conservation equation in the static case is (2-8b) from 
which the condition of hydrostatic equilibrium may be derived: 

^,rip + p)+P,r = 0. (3-4) 

The pressure is to be derived from this equation. The metric function A(r) is given 
directly by (2-10) 
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where it can be seen that the mass of the star increases as r D 1 . 

By using (2-11) and the equation of hydrostatic equilibrium, one obtains: 



a r + X(r) r — e X(r > = 0, 

' r w ' r 2 D — 2 



(3-6) 



with a := (p + po) 1 - 

After much calculation and consideration of boundary conditions at the 
surface of the star (the pressure must vanish at the stellar boundary), the above 
equation may be solved for the pressure: 



P = Po 



(D - 3)e" 



-A(r)/2 



(D - 3)e 



-A(a)/2 



(D - l)e- A ( a V 2 — (D — 3)e- A «/ 2 



(3-7) 



As described in the previous section, this is that last piece of information required 
to construct the solution. The metric function, v(r), is now computed directly from 
(2-12): 



- 9tt := = 



(D - l)e 



-A(o)/2 _ 



(D - 3)e 



-A(r)/2 



(D - l)e- A ( a )/ 2 - (D - 3) 
and total fluid mass of the sphere is simply 



M = 



8vrp 



(D-2){D-l) 



(3-8) 



(3-9) 



At this point it may be noted that the pressure becomes infinite at the same point 
at which the metric function, g u , vanishes. This happens when 



r 2 _{D- l) 2 Ma 2 - 2(D - 2)a 



D-l 



(3-10) 



M(L>-3) 2 

To ensure that this does not occur for any acceptable value of r, we must enforce: 



M < 2(D- 2) 



- (D- 1)2- 



(3-11) 



This is the upper bound on the gravitational potential of the Z)-dimensional star and 
is the generalisation of the M < |a law (Buchdahl's theorem 30 )) in four dimensions 
which applies to any spherical, static stellar model. 

Finally, the metric at the boundary takes on the form: 



ds \r=a 



=-A(a) 



2 dt 2 + e x(a) dr z +r z df2'( D _ 2) , (3-12) 



[{D - l)e~ x (a)/2 _ (£> _ 3 )] 
from which it can be seen that a coordinate re-scaling: 



(D - l)e 



-A(a)/2 _ 



(D-3) 



-l 



t, 



(3-13) 



allows smooth joining of (3T2) to the A = case of (3-2). 
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3.3. Anisotropic fluid 

Another example is that of the anisotropic fluid star characterized by: 



T^v = (p + P±)u^u u + p u g^ u + (p„ - p±)s fM s u , 
u a u a = — 1, s a s a = +1, tt Q s a = 0. 

In the static case, the metric (2-3) and (3-14) yield 

ds 2 = _ e „(r) ^2 + e A(r) ^2 + r 2 



(3-14) 



(3-15) 



n * = e -"(r)/2 r 



e -A(r)/2 S t = /(A) 



0. 



Now, the static versions of equations (2-10), (2-12) and (2-13) provide the general 
solution of the problem as: 



e -A(r) 



16vr 



(D - 2)r D ~ 3 
2m(r 



p(x)x u ~ 2 dx + 



-D-3 



(3.16a, 



''() 



„£>-3 



exp < c + 



16vr 



ro L 
L> - 2 



p„(x) +p(x) 
2m(x) 



D-3 



„-D-2 



3,- 16b) 



\v(r\ r + — 



=:p ± (3-16c) 



Here, fco and c are two arbitrary constants of integration. The constant c can be 
absorbed by the transformation i = e c l 2 t and the constant k can be set to zero to 
avoid a singularity at the center. 

3.4. D-dimensional Neutron star 

The neutron star represents a possible end state of a collapsed massive star. 
The electrons in the matter making up the star are subject to extreme pressures 
and, via the reaction 

p+ + e" -^n + u , (3-17) 

are converted to neutrons. The bulk of the energy is carried away by the neutrinos 
leaving behind a "cold" or degenerate remnant. 

The D-dimensional uncertainty principle per unit volume yields 



\[AV = h D -\ 

3=1 



(3-18) 



„th 



where k 3 is the momentum associated with the k spatial direction and h is the 
Planck's constant. By considering a spherical shell of inner radius k and thickness 
Ak in the phase space, the number of neutrons per unit volume in this shell will 
be the maximum occupation number times the number of cells with phase volume 
h D - 1 : 



4n 2 k 



D-2 



Ak 



(3-19) 
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By integrating (3-19) from k = to the Fermi momentum, kp, we get the maximum 
number of neutrons per unit volume with momentum up to /cp: 



AT 



47T V kp~ l 



r 



'D-v 



(D - l)^- 1 



The energy density of the system is simply given by (using (3-19)): 

-kp 



D-l 



r hp-* Jo 



k D ~ 2 (k 2 



m 



V 2 dk, 



(3-20) 



(3-21) 



which, in the extremely relativistic limit (acceptable for interior regions of neutron 
stars 15 )), yields: 



D-l n 



(3-22) 



The pressure of the system is given by differentiating the energy with respect 
to volume: 



dE 

P = ~dV 



D-l J-, 

4vr— k D 



djp/M) _ _ 

d{l/M) D{D- l)r (^i) h ' 1 D-l 



(3-23) 



This gives the equation of state. Using (3-23) in the conservation law (2-8b) and 
utilising (2-11) yields the equation of hydrostatic equilibrium for a D-dimensional 
neutron star (or any highly relativistic, degenerate D-dimensional Fermi gas): 



-D{D-3)p 



1 



2m(r) 

r D-3 



n -l 



L(D-3)(D-2)(D-1) 



+ m{r) 



(3-24) 



The densities in the above equation may be eliminated in favour of mass terms, 
giving an O.D.E. for the mass. By assuming a series solution for m(r), one finds 
that only one term in the series can contribute to the solution yielding: 



m(r) 



2(D - 1) 



,D-3 



D 2 (D-3) + 4(D-l) 



p(r) - (D-3)(D-2){D-1) r - 2 
^ ; 4vr[Z? 2 (Z)-3)+4( J D-l)] ' 



(3-25a) 
(3-25b) 



This second equation indicates that the fall-off properties of the density for the 
degenerate Fermi gas is independent of dimension. The metric functions may now 
be calculated via (2-10) and (2-12): 



u(r) 



1 



4{D - 1) 



D 2 (D - 3) + 4(D - 1)' 
4(D - 1) 



D 2 (D-3)+4(D-l) 



{r/r )D . 



(3-26a) 
(3-26b) 



Note that for all D > 3, the metric given by (3-26a) and (3-26b) is an i?-domain 
type metric. 
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The above metric derivation is only valid for inner layers of neutron stars 
where one can make the extreme relativistic approximation. The singularity in the 
density at the center of the star is also a manifestation due to this approximation. 
This singularity plagues the known four dimensional model as well where the region 
near r = is usually excised. 

§4. The existence of Birkhoff's theorems 

One very interesting aspect of spherical symmetry is the existence of Birkhoff's 
theorem 31 ) . The original theorem applies to four dimensional systems and states that 
spherically symmetric vacua are static and are locally equivalent to the Schwarzschild 
solution. This theorem has since been generalised to include certain matter fields in 
4 dimensions including electromagnetic 32 -* , 33 ^ and scalar fields 34 ) , 35 ) . A study of the 
differentiability properties required for a well posed staticity theorem may be found 
in 36 ). Also, some elegant techniques have been employed in the literature regard- 
ing higher dimensional versions of the theorem 37 ) , 38 \ We briefly present here the 
conditions required for, as well as a proof of, a D-dimensional Birkhoff's theorem 
which applies to both vacuum and non-vacuum systems. The theorem presented 
here generalises previous theorems and below we will extend the theorem to cases 
where the metric is C° (piece-wise C°). The assumptions required for the general 
theorem to hold are minimal and, most likely, cannot be relaxed, allowing for a most 
general theorem. We do, however, examine specific cases where conditions may be 
weakened and the theorem will still hold. We state the theorem as follows: 

Theorem: 1 Let B C R 2 be a convex domain in the r — t plane. Let spherically 
symmetric metric functions g rr > and gu < belong to the class C 3 (B) and the 
stress-energy tensor, T^, belong to the class C 1 (B). Moreover, let: i) T l r = and 
ii) T r rt = 0. Then, the metric solutions satisfying (2-7a)-(2-7d) must admit an 
additional Killing vector. 

Proof: By the assumption on the metric functions we can write gu(r,t) = —e u ( r '^ 
and g rr (r,t) = e A ( r '*) where f{r,t) and A(r, t) are of class C 3 . The identity T\, = 
yields, from the equation (2-7c), that 

A = A(r). (4-1) 
From this, the equation (2-7a) implies that 

T\ = T\{r). (4-2) 

Utilising (4-1) in (2-7b) yields: 

Kr , () , _ ^ {s, nW + < D - 2 >'f- 3 » [l - .-*>] } . ,4-3) 
Therefore, by differentiability with respect to t, we obtain 



v(r,t), r ,t = 



(4-4) 
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(D-2)- 



(4-5) 



(4-6) 



By a coordinate transformation, 




(4-7) 



the metric in (4-6) reduces to a static one which admits the additional Killing vector 



Note that the only physical assumptions involved in the theorem are staticity of 
the radial pressure and that the system possess no mechanism for radial energy 
transport. In some circumstances, the above assumptions may be weakened as will 
be discussed below. 

Convexity of the domain is required to address certain situations where the 
theorem fails though it should seemingly otherwise hold. Consider the following 
simple counter-example 39 ^ depicted in figure 1. 



Fig. 1. The non-convex domain provided by considering the area bounded by the the dashed 
line. The line L is removed from the region creating a non-convex domain. In this case, 
a staticity theorem will not hold. 

The line segment, L, is given by L := {(r, t) G M. 2 : 1 < r < 2; t = 0}. The 
non-convex domain is furnished by considering the bounded rectangular domain with 
the line, L, removed: B := (0,2) x (—2,2) — L. Consider now the metric function, 
X(r, t), given by 



d_ 



(3/2, 3/2) 



1 




(3/2, -3/2) 



J 




for 

otherwise 



2 > r > 1, t > 
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a C function with derivatives: 

X(rt) -{ (r " 1)4 f ° r 2>r>l, t>0 . . 

Air ' tj ' r " \ otherwise 14 9 j 

and 

\(r,t), t = 0. (4-10) 

Consider the two points (|, |) and (|, ^) in the r — t plane. The function 
A(r, t) possesses values: 



A 



3 3\ _ 1 
2' 2 J ~ 160 ' 



*(H) =» 

Therefore, although A(r, t) $ = in the domain, A(r, i) is no£ independent of time and 
the metric is non-static. Such arguments may be applied to cases where non-trivial 
topological features in the manifold will create non-convexity in the r — t domain. 

The conditions under which the original, four dimensional vacuum, theorem 
hold have been weakened to the point of admitting a C° metric 36 ) if the metric 
possesses a separable gu- We weaken further the conditions here and demonstrate 
that a D-dimensional metric (not necessarily vacuum) with separable gu may possess 
behaviour as pathological as piece-wise C° (a gravitational shock front) and still be 
static. 

With little loss of generality, the metric under consideration may, in the 
vicinity of the "jump", be written as 

ds 2 = _ e a(r) B 2^ ^2 + g A(r) ^2 + r 2 df2 f D _ 2) , (4-12) 

where 

B{t) := 2 + Sgn(t) . (4-13) 

The function Sgn(t) is given by 

%■*>=-{ T t Zl ^ 

One may believe that such a shock front would create singular structure in the 
manifold since 

dB(t) 
dt 

as well as 



2 <$(*), (4-15) 



r* = Un{-g tt ) t . 

rit = \a(r)^-^B\t). 

Surprisingly, no singularity is present in the curvature tensor and such a solution 
possesses well behaved orthonormal Einstein tensor. This is because all derivatives 
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of g tt are multiplied by quantities which vanish at t = and it has been rigorously 
proved that • 5(t) = 40 \ The transformation in (4-7) is equally well defined as: 

t = f B{r)dT t < 
J t 

t + \t Q \ for i < 



3t+\t \ for t>0. ^ 4 ' 16 ) 

Notice that the re-scaled time variable, t, is continuous at at t = 0. Such a trans- 
formation is therefore admissible and the re-scaled metric is explicitly static. It is 
unknown if behaviour worse than C° (such as characteristic fluctuations) may admit 
a staticity theorem. 

It is well known that hyperbolic equations admit exact discontinuous solu- 
tions and that Einstein's equations therefore allow for such solutions. Before t = to 
the metric under consideration can be transformed smoothly to a static one. After 
t = to a similar transformation also renders a static metric. At t = to there exists 
a C° transformation to a static metric. Therefore, in the entire domain, the metric 
can be transformed to a static one although the transformation is not the usual C 3 . 

The discontinuous metric in the form of (4T2) yields discontinuous or- 
thonormal tetrad as well as the singular Christoffel symbols mentioned above. How- 
ever, the Riemann invariants, which govern the gravitational physics possess a re- 
movable dicontinuity which may be ignored. The potentially problematic compo- 
nents are: 



2 



"M,r,r + ^a 2 r - -a(r) )r .A(r), 



where hatted indices denote quantities calculated in the orthonormal frame. Notice 
that the Riemann invariants are well behaved and independent of the function B(t). 
These components are all static indicating staticity of the gravitational field and 
therefore yielding a Birkhoff's theorem. 

Generally, a Cq metric yields serious singularities in the manifold (see 43 -*, 
44 ) for studies of such systems). However, for the type of metric given by (4-12), 
singularities are avoided. 

Finally, we should mention that results in this section would be unaltered if 
one relaxes the condition for the space-time to possess a strictly Lorentzian metric. 
In such case, the metric (4-6) is better written as 

ds 2 = - e «Mi(t) dt 2 + e A W dr 2 + r 2 df2f D _ 2) , (4-17) 

where the function A(t) may switch sign at one or more values of t, such as A(t) = t 3 
(see 41 -* 42 ) for examples of such space-times). The Birkhoff's theorem holds both in 
the Lorentzian branch and the Euclidean branch of the manifold. Also, in the case 
where the signs of gtt and g rr are both switched, we obtain the generalisation of 
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Birkhoff's theorem in the T-domain 45 ) . In this situation, the theorem is a homo- 
geneity theorem as opposed to a staticity theorem since the additional killing vector 
is space-like. 

§5. Concluding remarks 

We have considered spherically symmetric solutions to the Einstein field 
equations with an arbitrary number of spatial dimensions. A reasonably general 
method has been presented which allows one to solve, at least in quadrature, these 
equations. The method has been illustrated by quickly and efficiently computing 
the metric for a D-dimensional black hole with arbitrary cosmological constant, 
the metric for a D-dimensional constant density star, anisotropic fluid star and 
neutron star. An upper limit has been placed on the mass/radius ratio of stars of 
arbitrary dimension. Finally, the minimum general requirements for a D-dimensional 
Birkhoff's theorem in both vacuum and non-vacuum systems has been presented. To 
have a rigorous theorem, one must insist on convexity of the domain in question. In 
certain situations, a staticity theorem may hold even when the metric component, 
gu, is only C® or piece-wise continuous. It may be shown that the theorem also holds 
for metrics of Euclidean signature as well as inside black hole T-domains. 
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